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Numerical and semi- analytical models are presented for photon-enhanced-thermionic-emission (PETE) de- 
vices. The models take diffusion of electrons, inhomogeneous photogeneration, and bulk and surface recom- 
bination into account. The efficiencies of PETE devices with silicon cathodes are calculated. Our model 
predicts significantly different electron affinity and temperature dependence for the device than the earlier 
model based on a rate-equation description of the cathode. We show that surface recombination can reduce 
the efficiency below 10 % at the cathode temperature of 800 K and the concentration of 1000 suns, but 
operating the device at high injection levels can increase the efficiency to 15 %. 



I. INTRODUCTION 

Currently solar energy is converted to electric power 
using two technologies: photovoltaic (PV) solar cells and 
concentrated solar power systems based on heat engines 1 -. 
The former system requires low and the latter high oper- 
ating temperatures. This discrepancy poses a challenge 
for combination of the two systems in tandem, where the 
heat engines exploit the waste heat of the PV system. 
A photon-enhanced-thermionic-emission (PETE) device 
proposed by Schwede et aL - is a PV device which bene- 
fits from high operation temperatures. It can be coupled 
to a heat engine, thereby allowing total efficiencies above 
50% to be potentially reached. 

The PETE device is depicted in Fig. Q] The photons 
are absorbed in the cathode, i.e. the absorber, which is a 
P-type semiconductor. The cathode material should have 
a suitably low energy gap so that most of the solar pho- 
tons are absorbed. The absorbed photons are themion- 
ically emitted from the cathode to vacuum, where they 
travel to the anode, i.e. the electron collector. The sur- 
face of the cathode should have a low electron affinity 
in order to have reasonably strong thermionic emission. 
Electron affinity can be tuned significantly below the bulk 
value by different surface coatings (see Refs. [H and [IH 
and references therein). Even negative electron affinities 
can be obtained for silicon 1 ^. These coatings, however, 
might not be stable at the high operation temperatures 
of PETE devices. The anode material can be metal or 
N-type semiconductor with suitably low work function in 
order to have high output voltage V for the device. 

In the PETE device model of Ref. the cathode ma- 
terial is described by a single rate equation that neglects 
many important effects, such as diffusion and realistic 
recombination. In this article, we present a model that 
takes the most of the relevant effects in semiconductors 
into account. We solve the electron density in the cath- 
ode numerically in the general case and derive also a 



semi-analytical model, which can be used at low injec- 
tion levels. Silicon is a good candidate for PETE due 
to its rather low band gap, good thermal stability, high 
availability, cost-efficiency, and good manufacturability. 
Therefore, we use our models to calculate the character- 
istics of a PETE device with a silicon cathode using ex- 
perimentally verified material data. We find that at low 
temperatures our model predicts a significantly lower ef- 
ficiency than the model presented in Ref. |2|. The overall 
temperature and electron affinity dependency of the effi- 
ciency differs also from that of Ref. |2|. Furthermore, we 
show that surface recombination can reduce the efficiency 
of the PETE device below 10 %, but operating the de- 
vice at high injection levels can provide an enhancement 
where efficiency of 15 % is reached. 



II. THEORY 

A. Semiconductor material model 

The Fermi level Ep is solved numerically using the 
electroneutrality condition- n eq + = p eq + , where 



'eq 



and p eq are the densities of electrons and holes in the 
thermodynamical equilibrium, and and Np are the 
densities of ionized acceptors and donors, respectively. 
For the temperature dependence of the band gap of sili- 
con we use the standard formula 3 



E g (T) = E g (0) 



T + f3 g > 



(1) 



where E g (0) = 1.170 eV, a g = 4.73 x 10" 4 eV/K, and 
f3 g = 636 K. 

The total minority-electron lifetime in bulk can be 
written as 4 ^— 



1 



B(n+p eq ) + (C n n + C p p)(n + p eq ) + 



(2) 
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where B is the radiative recombination coefficient, n 
and p are the densities of electrons and holes, respec- 
tively, C n and C p are the Auger recombination coef- 
ficients for electrons and holes, respectively, and tsrh 
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FIG. 1. (a) Schematic picture of the PETE device. T c is 
the absolute temperature of cathode, Ta the absolute tem- 
perature of anode, V the output voltage, J the density of 
the output current, and R a resistor representing the external 
load. The energy bands of the device (b) in the thermody- 
namical equilibrium and (b) in the flat-band case, which often 
yields the highest power output. Ef is the Fermi level, Ef,v, is 
the quasi Fermi level of electrons, q is the elementary charge, 
Xc is the electron affinity of the cathode, <f>c and 4>a are the 
work functions of the cathode and the anode, E v is the va- 
lence band maximum, E c is the conduction band minimum, 
and E vac is the vacuum energy level. Here the anode is metal, 
but it can be an N-type semiconductor as well. 



is the Shockley-Read-Hall (SRH) recombination lifetime, 
which is directly proportional to the density of the deep- 
level impurities in the material, and it increases when 
the injection level is increased 5 -. The SRH lifetime has 
constant saturation values both in low and high injec- 
tion regimes 5 . Therefore, we use constant tsrh = 2.5 /is 
to represent a high quality silicon^ at either low or high 
injection level. For silicon 5 - B = 4.73 x 10~ 15 cm 3 /s and 
C p = 10 ~ 31 cm 6 /s, which is a value measured at low 
injection level. In a P-type material C p determines the 
Auger lifetime in low and the ambipolar Auger coefficient 
Ca = C n + Cp in high injection conditions 5 , respectively. 
The experimental high-injection value for Ca differs from 
C n + C p measured at low injection levels^. We take 
the both injection regimes into account using an artifi- 
cial value of 10~ 30 cm 6 /s for C n . For electron mobility fi n 
in silicon we use a modelZ- optimized for a wide tempera- 
ture range. Mobility is linked to the diffusion coefficient 
D n = ksT Hn/ q, where T is the absolute temperature 
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FIG. 2. Schematic picture of the cathode of the PETE device. 



and q the elementary charge. For the absorption coef- 
ficient a of silicon we use a widely-used semi-empirical 
models. 



B. Electron density 

The electron density n in the cathode sketched in Fig. [2] 
can be calculated using the continuity equation, the gen- 
eration and recombination rate equations, and the drift- 
diffusion current formulas for the charge carriers 3 -. We 
assume that the electric field inside the cathode is ap- 
proximately zero, because the most of the voltage differ- 
ence is across the vacuum gap and the photogeneration of 
electron-hole pairs supports the local approximate charge 
neutrality. This simplifies the solution considerably, since 
holes need not to be taken into account explicitly. In ad- 
dition, the band bending near the cathode surface (see 
Fig. [Tp) is neglected, since the effect is small and the 
PETE device is mostly operated near the flat-band case 
(shown in Fig. QJ), where the band bending disappears. 
The excess-electron density An(x) — n(x) — n eq can be 
described by 



d 2 An 
dx 2 



An _ G n (x)T„ 
1J 



LI 



(3) 



where L n = y ' D n T n is the diffusion length of electrons 
and the generation is given by 



G(x) = a$ (1 - R FS ) 



+ RBse 



a(x-2W) 



(4) 



where a is the absorption coefficient, $n is the incident 
flux of photons, Rfs(bs) is t ne front (back) surface re- 
flection coefficient, and W is the thickness of the cathode. 
The boundary conditions are 
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where 5' 



nFS(BS) is the front (back) surface recombination 
velocity and J is the density of the output current. 

The analytical solution exist when r„ is constant, 
which is valid at low injection levels (i.e. An Peq)- 
Assuming the low-injection case and using Eqs. ([5]) and 
^ the solution of Eq. at x — W is given by 
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where A is the wavelength of photons. The first term 
in Eq. ([7]) describes the fact that the electrons need to 
diffuse from bulk of the cathode to the back surface in 
order to be emitted. When current is drawn, An(W) will 
decrease because of this diffusion process. The second 
term in Eq. ([7]) represents the excess electrons due to the 
photogeneration. 



C. Output current 

The density of the cathode current, which corresponds 
to electrons emitted from the cathode, can be written 
using the quasi Fermi level Ep >n = Ep + k B Tc ln (n/n eq ) 
as 



J c = A£T£exp 



AE C \ n 
ksTc I n eq 



(12) 



where AEc — q4>C when V < Vfb, and AEc = q&c + 
q(V — Vfb), when V > Vfb, A* c is Richardson's constant, 
and Tc is the absolute temperature of the cathode. The 
flat-band voltage Vfb is defined as 



Vfb = 4>c - 4>a, 



(13) 



where 6 a is the work function of the anode material and 



t>c=-(E c -E F +xc) 

q 



(14) 



is the work function and xc the electron affinity of the 
cathode material. The density of the anode current, 
which corresponds to electrons emitted from the anode, 
is given by 



J A = A* A T%exp [ - 



AEa_\ 
ksTA J ' 



(15) 



where AEa — q4>A + q{Vfb — V) when V < Vfb, and 
AEa — q4>A, when V > Vfb, A* A is Richardson's constant 
of the anode, and Ta is the absolute temperature of the 
anode. We assume that all the electrons emitted from 
the cathode are collected by the anode and vice versa: 



J = Jc — J A- 



(16) 



The flat band voltage Vfb is an important parame- 
ter for the efficiency of the PETE device. At voltages 

V above Vfb the additional energy barrier q(V — Vfb) 
for the electrons emitting from the cathode appears (see 
Eq. (fT2"jn . This decreases Jc considerably. At voltages 

V below Vfb the electrons emitted from the anode are 
hindered by the energy barrier q(Vfb — V) (see Eq. |[T5])). 
This allows J a to be reduced by lowering V . In gen- 
eral, the highest output power will often be obtained near 

V = Vfb- However, when Tc ^> Ta, the high thermal en- 
ergy of the cathode electrons allows the range V > Vfb 
to be used as well. 

Using n eq = 7V c exp [-(E c - E F )/(k B T c )} and Eqs. © 
and P^)) - (fTB|) the output electric current density can be 
written alternatively as 



J '^Sl 



ark i 



where the photocurrent is given by 
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where 
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for V < Vfb 
k for V > V fb . 



The dark current is given by 



^dark — 



qD n n, 



eq 



dL r . 



A\T\ {qAE dark 



A* c Tl 



(19) 



(20) 



where A^ dark = 6 C + (V - 6 C )T C /T A when V < V fb , 
and A£; dark = V + 6 A {1 - T C /T A ), when V > V fb . The 
dark current is the electric current that flows through 
the device when there is no illumination. Under illu- 
mination it usually decreases the total current, thus it 
should be eliminated. This can be done by choosing V 
optimally. However, if Ta is very low, the direction of 
■/dark can change, and the device harvests energy also 
from the thermal energy of cathode electrons similarly 
as a thermionic converter. 

Eqs. (fT8|l and (fT9|) show that in order to maximize 
J sun xc an d k should be as small. On the contrary, then 
>/dark will also increase. J sun can be increased and Jdark 
decreased by increasing Tc- Jdark can be reversed by 
having high 6 a and low &c (this will reduce Vfb)- In 
addition, T A and A* A should be small and Tc and A* c 
should be large. 
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FIG. 3. Efficiencies of PETE devices with silicon cathodes 
with various electron affinities X c as functions of the cathode 
temperature Tc calculated using the numerical model with 
Na = 10 18 cm -3 and W = 5 j«m and without surface recombi- 
nation at the concentration of 1000 suns. The results from the 
rate-equation model 2 - were calculated with E g (T) = 1.12 eV. 
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FIG. 4. Efficiency of a PETE device with a silicon cath- 
ode and various back-surface recombination velocities S u bs as 
function of the cathode temperature Tc calculated using the 
numerical (solid lines) and semi-analytical (dotted lines) mod- 
els with SnFS = 10 2 cm/s, X c = 0.4 eV, N A = 10 18 cm -3 , 
and W = 5 fim at the concentration of 1000 suns. 



III. RESULTS AND DISCUSSION 

The semi-analytical model, Eqs. ©, pl]l. and (fT71) - 
(I2TJ1) . applies in the low injection condition (An <C p eq ) 
with a P-type cathode material (jp eq ^> n eq )-. The nu- 
merical calculations were performed using Eqs. ([3])-© 
and (IT2|) - (|T6)) . In both models V was optimized nu- 
merically for maximum output power. The AMI. 5 di- 
rect+circumsolar spectrum applicable to solar concentra- 
tors was used in the calculations. We use (111) silicon as 
the cathode surface, for which 3 A* c — 264 A/cm 2 /K 2 . 
For the material parameters of the anode we use the 
same values as Schwede et al.— , 4>a = 0.9 V and A* A = 
120 A/cm 2 /K 2 . T A was set at 573.15 K in order the heat 
engine potentially coupled to the anode to have a reason- 
ably high efficienc y 10 ' 11 . We also assume Rfs = and 
Rbs = 1 for simplicity. The effect of Rbs, however, is 
not very large (see below), since it has an effect only on 
the absorption of low energy photons. 

The efficiencies of PETE devices with various values of 
Xc are plotted as functions of Tc in Fig. [3] They increase 
with increasing Tc and decreasing xc due to the increase 
of Jc- Using Rbs = instead reduces the efficiency from 
13.0 % to 11.7 % at 550 K in the case xc = 0.4. 

Fig. [3] shows also the efficiencies given by the simple 
rate-equation model 2 -, which assumes that all photons 
with energy greater than E g are uniformly absorbed in 
the cathode. Only the uniform radiative recombination 
is taken into account with a general model based on the 
black-body radiation. At low Tc the rate-equation model 
suggests much higher efficiencies than the present more 
complete model. This is mostly due to the facts that all 
the possible photons are absorbed and the Auger recom- 
bination is not included in the rate-equation model. 

The efficiency is a balance between many opposing ef- 
fects, which depend on Na and W: Increasing Na in- 
creases Vfb which increases the efficiency. On the other 
hand, high Na reduces the efficiency due to the Auger 



recombination which is proportional to N\. The effect 
of the bulk recombination can be reduced by decreasing 
W, but then some of the photons will not be absorbed. 
Although 5 /im is already small thickness for a silicon 
absorber, but, if W is increased to 50 /im while keeping 
Na constant, the efficiency reduces from 13.0 % to 8.4 % 
at 550 K in the case xc = 0.4 eV. At high T c the ef- 
ficiency curves of the present model unite regardless of 
the differences in xc- The reason for this is that the to- 
tal current does not depend on xc in this range: When 
Tc ~> Ta the range V > Vfb can be utilized and the out- 
put current density can be written using Eqs. (|12[) - (fT(3|) 
as J = A c T c exp[(E F , n - E F - qV - qf A )/(k B T c )\ - 
A* A T A exp(—q(/)A/kBTA)- The rate-equation model does 
not behave like this at high Tc, because An decreases 
much faster with increasing Tc than in the present model. 
This is caused mostly by the differences in the model- 
ing of recombination and the narrowing of the band gap 
with increasing Tc which is taken into account only in 
the present model. The latter effect enhances both pho- 
togeneration and thermal generation of electrons, which 
causes the apparent efficiency increase at Tc > 800 K. 

The effect of S n BS on the efficiency of the PETE de- 
vice is shown in Fig. 0) Even a low value of S n ss re- 
duces the efficiency remarkably. Similar results were also 
obtained with various values of S u fs- The efficiency is 
more sensitive to S u fs than S n BS because photogener- 
ation is much stronger near the front surface than near 
the back surface. The semi-analytical model is in a very 
good agreement with the numerical model when the re- 
combination velocities are > 1000 cm/s since the injec- 
tion level is below 0.1 at Tc > 550 K. In the low-injection 
conditions achieving surface recombination velocities less 
than 100 cm/s requires usually use of back surface field 
structures 3 -. However, the surface recombination velocity 
decreases rapidly when the injection level is increased^ 
and values below 2 cm/s can easily be reached at injec- 
tion level of unity (details depend on the properties of 
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the surface states} 5 . Therefore, the case of zero surface 
recombination can actually be realistic for a PETE de- 
vice. 

IV. CONCLUSIONS 

In summary, we have built a theoretical model for 
PETE devices. Our model takes electron diffusion, inho- 
mogeneous photogeneration, and bulk and surface recom- 
bination into account. In comparison to the rate equa- 
tion model of Ref. our model predicts different depen- 
dency of efficiency on such paramaters as cathode elec- 
tron affinity and temperature. In most cases our model 
also predicts lower efficiency. Especially, the surface re- 
combination present on real surfaces can reduce the ef- 
ficiency to extremely low values. However, the surface 
recombination might have only a very weak effect on the 
performance, since the PETE device works often in the 
high-injection conditions where the surface recombina- 
tion can be rather small 5 -. We finally point out that the 
realization of the PETE device requires choices of many 
parameter values and materials which should all be opti- 
mized. In addition, the heat balance between the PETE 
device and the heat engine coupled to the anode, which 
was not considered in this article, should be managed as 
well. Full optimization with heat balance modelling will 
be left for future studies. 
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